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SUMMARY 


This  analysis  treats  the  air  forces  and  moments  in  supersonic 
potential  flow  on  oscillating  triangular  wings  and  a  series  of  swept- 
hack  and  arrow  wings  with  subsonic  leading  edges  and  supersonic  trailing 
edges.  For  the  wings  undergoing  sinusoidal  torsional  oscillations  simul¬ 
taneously  with  vertical  translations,  the  linearized  velocity  potential 
is  derived  in  the  form  of  a  power  series  in  terms  of  a  frequency  param¬ 
eter.  This  method  can  be  useful  for  treatment  of  similar  problems  for 
other  plan  forms  and  for  wings  undergoing  other  sinusoidal  motions. 

For  triangular  wings,  as  many  terms  of  such  a  series  expansion  as  may 
be  desired  can  be  determined;  however,  the  terms  after  the  first  few 
become  very  cumbersome. 

Closed  expressions  that  include  the  reduced  frequency  to  the  third 
power,  an  order  which  is  sufficient  for  a  large  class  of  practical 
applications,  are  given  for  the  velocity  potential  and  for  the  components 
of  chordwise  section  force  and  moment  coefficients. 

These  wings  are  found  to  exhibit  the  possibility  of  undamped  tor¬ 
sional  oscillations  for  certain  ranges  of  Mach  number  and  locations  of 
the  axis  of  rotation.  The  ranges  of  these  parameters  are  delineated  for 
triangular  wings. 


INTRODUCTION 


This  paper  is  concerned  with  the  derivation  of  expressions  for  the 
velocity  potential  and  associated  forces  and  moments  for  oscillating 
triangular  wings  in  supersonic  flow.  The  boundary-value  problem  for 
the  linearized  velocity  potential  for  an  apex-forward  triangular  wing 
oscillating  in  a  supersonic  main  stream  may  be  classified,  according  to 
reference  1,  as  "purely  supersonic"  if  the  leading  edges  of  the  triangle 


are  outside  the  Mach  cone  emanating  from  the  apex  of  the  triangle  or 
T,mixed  supersonic"  if  the  leading  edges  are  inside  this  Mach  cone. 

In  the  purely  supersonic  case  the  principle  of  independence  holds ; 
that  is,  the  flow  on  the  upper  surface  of  the  wing  is  independent  of  the 
flow  on  the  lower  surface  and  vice  versa.  Garrick  and  Rubinow  (refer¬ 
ence  l)  have  shown  that  the  boundary -value  problem  for  the  velocity 
potential  in  the  purely  supersonic  case  can  be  satisfied  by  simple  dis¬ 
tributions  of  sources  with  local  strength  proportional  to  the  local 
prescribed  normal  velocity  of  the  wing. 

In  the  mixed  supersonic  case  the  principle  of  independence  does  not 
hold.  Boundary -value  problems  for  lift-producing  wings  in  this  case  can 
be  satisfied  by  distributions  of  doublets;  the  relation  between  doublet 
strength  and  normal  velocity  of  the  wing  is,  however,  in  general,  not 
simple.  The  determination  of  this  relation  requires  the  solution  of  an 
integral  equation  that  employs  the  potential  of  a  time-dependent  unit 
doublet  as  kernel  and  limits  of  integration  that  depend  on  Mach  number 
and  wing  plan  form. 

For  treatment  of  problems  that  involve  boundary  conditions  that  are 
independent  of  time,  such  as  constant  angle  of  attack,  constant  rate  of 
pitching,  and  so  forth,  the  doublet  potential,  employed  as  kernel  of  the 
integral  equation,  is  considered  independent  of  time  and  in  these  cases 
the  integral  equations  for  triangular  wings  can  be  solved  by  a  straight¬ 
forward  process. 

For  treatment  of  problems  of  oscillating  wings,  however,  it  is 
necessary  to  employ,  as  the  kernel  of  the  integral  equations,  a  doublet 
potential  that  varies  harmonically  with  time  and  in  this  case  the  solu¬ 
tion  of  the  integral  equation,  generally,  becomes  very  cumbersome.  If 
the  doublet  potential  or  kernel  is  expanded  in  terms  of  the  frequency 
of  oscillation,  however,  use  can  be  made  of  knowledge  of  solutions  of 
integral  equations  for  problems  that  are  independent  of  time  to  obtain 
an  expanded  form  of  solution  for  a  wing  undergoing  harmonic  oscillations. 
Such  a  procedure  was  demonstrated  in  treatments  of  rectangular  wings  in 
references  2  and  3 •  (The  derivation  in  reference  3  is  based  on  an 
erroneous  argument,  regarding  certain  terms  in  the  normal  velocity; 
nevertheless,  the  final  expression  given  for  the  velocity  potential  is 
correct. ) 

The  purpose  of  the  present  paper  is  to  make  use  of  the  expanded 
form  of  the  velocity  potential  to  obtain  the  forces  and  moments,  based 
on  the  first  few  terms  of  this  potential,  for  a  rigid  triangular  wing 
performing  vertical  and  pitching  sinusoidal  oscillations  in  mixed  super¬ 
sonic  flow.  Although  as  many  terms  of  the  expanded  potential  as  may  be 
desired  can  be  obtained  after  the  first  few  terms,  the  process  becomes 
very  cumbersome.  The  flow  normal  to  the  leading  edge  is  subsonic  but 
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the  flow  normal  to  the  trailing  edge  is  considered  to  he  supersonic. 

This  latter  consideration  implies  that  the  potential  derived  for  trian¬ 
gular  plan  forms  may  be  used  to  calculate  the  aerodynamic  forces  and 
moments  for  other  plan  forms  that  may  be  formed  with  the  triangular  wing 
by  cutting  the  trailing  edges  so  that  they  lie  ahead  of  the  Mach  cones 
emanating  from  their  foremost  points. 

Other  approaches  to  the  solution  of  the  problem  of  oscillating 
triangular  wings  have  been  given  by  Robinson  (reference  4),  Haskind  and 
Falkovich  (reference  5)>  and  by  Stewartson  (reference  6).  In  both 
references  4  and  5  formal  solutions  to  the  problem  were  obtained  in  terms 
of  special  systems  of  curvilinear  coordinates.  Robinson's  solution  was 
given  in  terms  of  a  double  summation  of  trilinear  combinations  of  Bessel 
functions  of  the  first  kind  with  Lame'  functions  of  the  first  and  second 
kinds.  Similarly,  the  solution  of  Haskind  and  Falkovich  was  given  in 
terms  of  summations  of  Bessel  functions  of  the  first  kind  combined  with 
elliptic  integrals  of  the  first  and  second  kinds.  In  both  references  4 
and  5  the  potentials  were  not  reduced  to  useful  forms  for  calculating 
forces  and  moments. 

In  reference  6  Stewartson  makes  an  interesting  though  specialized 
use  of  the  Laplace  transformation  to  develop  a  method  whereby  terms  of 
the  velocity  potential  for  triangular  plan  forms,  expanded  as  herein, 
can  be  obtained.  Stewartson  gives  formulas  that,  except  for  errors 
presumably  in  printing,  can  be  used  to  develop  the  potential  to  the 
second  power  of  the  frequency  but  he  omits  many  details  in  his  derivation. 


SYMBOLS 


0  disturbance-velocity  potential 

x,y,z  rectangular  coordinates  attached  to  wing  moving  in 

negative  x-direction 

|,tj  rectangular  coordinates  used  to  represent  space  location 

of  doublets  in  xy-plane 

function  defining  mean  ordinates  of  any  chordwise  section 
of  wing  such  as  y  =  yp  as  shown  in  figure  1 

w(x,yi,t)  vertical  velocity  at  surface  of  wing  along  chordwise 
section  at  y  =  yp 

*6 


abscissa  of  axis  of  rotation  of  wing  as  shown  in  figure  1 
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t 

time 

k 

h 

vertical  displacement  of  axis  of  rotation 

4? 

*0 

amplitude  of  vertical  displacement  of  axis  of  rotation, 
positive  downward 

a 

angle  of  attack 

“o 

amplitude  of  angular  displacement  about  axis  of  rotation, 
positive  leading  edge  up 

•  * 
h,a 

time  derivatives  of  h  and  a,  respectively 

V 

velocity  of  main  stream 

c 

velocity  of  sound 

M 

free-stream  Mach  number  (v/c) 

3  =  \jyi2  -  l 

* 

0) 

frequency  of  oscillations 

5,i% 

vp2 

k 

reduced  frequency 

e 

half  apex  angle 

C  =  tan  e 

anm 

represents  functions  of  oo,  x,  and  M 

functions  used  to  denote  doublet  distribution  functions 

dn 

constants  associated  with  Dn  depending  on  PC 

Ai'*\) 

constants  depending  on  PC 

°k 

constants  depending  on  PC  and  M 

m' 

2b 

root  chord  of  wing 

m 
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k 

* 


Q 

a 

q,K,6yr 

P 

Ap 

P 


semi span  of  wing 

slope  of  ray  passing  through  vertex  of  wing 

dummy  variables 

density 

local  pressure  difference 

section  force  (total  force  at  any  spanwise  station) 
components  of  section  force  coefficients 


Ma  section  moment  (total  moment  about  x  =  xq  at  any 

spanwise  station) 

Mi,M2,M3,M4  components  of  section  moment  coefficients 
%  total  component  of  damping-moment  coefficient 


ANALYSIS 

Boundary -Value  Problem  for  the  Velocity  Potential 


Referred  to  a  rectangular  coordinate  system  moving  forward  at  a 
uniform  supersonic  speed  in  the  negative  x-direction  (see  fig.  1)  the 
differential  equation  for  the  propagation  of  small  disturbances  that 
must  be  satisfied  by  the  velocity  potential  is 


C2\Bt  dx 


dx2  dy2  dz2 


(1) 


The  main  governing  boundary  condition  to  be  satisfied  by  the  velocity 
potential  is  that  the  flow  be  tangent  to  the  surface  of  the  wing,  or 


=  w(x,y,t) 


SZm  c)Zm 

dx  St 


where  2^  is  the  vertical  displacement  of  any  point  of  the  wing.  For 
the  particular  case  of  a  wing  independently  performing  small  sinusoidal 


6 


NACA  TN  2457 


pitching  oscillations  of  amplitude  Oq  about  some  spanwise  axis  xq 
and  small  sinusoidal  vertical  translations  of  amplitude  hQ,  the 
quantity  in  equation  (2)  is 

Zm  =  elcut|a0(x  -  xq)  +  h^J  =  a(x  -  xq)  +  h  (3) 

(See  fig.  l(b)  for  sketch  showing  instantaneous  displacement  of  sec¬ 
tion  y  =  y]_.)  For  convenience,  the  frequency  of  oscillation  of  both 
pitching  and  translatory  motion  is  denoted  by  a).  Considering  these 
motions  to  occur  at  separate  frequencies  would  add  no  difficulties  to 
the  derivation. 

Substituting  the  expression  for  2^  (equation  (3))  into  equation  (2) 

gives 


•  • 

w(x,y,t)  =  Va  +  a(x  -  xq)  +  h  (4) 

Equation  (4)  implies  that  the  velocity  potential  may  be  expressed  as  the 
sum  of  separate  effects  due  to  position  and  motion  of  the  wing  associated 
with  individual  terms  of  this  equation,  namely 

0  =  0a  +  0a  +  0h  (5) 


Derivation  of  0 

In  order  to  obtain  the  analytical  expression  for  the  potential  0, 
it  is  necessary  to  derive  only  one  of  the  subsidiary  potentials  appearing 
in  equation  (5),  say  0a.  The  other  subsidiary  potentials  0^  and  0£ 

can  then  be  obtained  from  the  derived  expression  for  0a  by  simple 
comparison. 

In  order  to  satisfy  the  boundary- value  problem  for  0a,  a  convenient 

procedure  is  to  start  with  the  expanded  form  of  the  potential  of  a  uni¬ 
form  distribution  of  doublets.  Then,  for  a  given  power  of  the  frequency 
of  oscillation  this  potential,  as  will  be  shown  in  the  following  analysis, 
can  be  modified  so  that,  when  integration  is  made  over  the  appropriate 
region,  the  results  satisfy  the  differential  equation  (l)  to  the  given 
power  of  the  frequency  and  satisfy  the  condition  of  tangential  flow 
exactly.  The  type  of  doublet  required  is  that  with  its  axis  normal  to 
the  plane  of  the  wing.  The  potential  of  such  a* doublet  may  be  obtained 
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from  the  potential  of  a  source,  located  in  the  plane  of  the  wing,  by 
partial  differentiation  with  respect  to  the  direction  normal  to  the 
plane  of  the  wing.  Similarly,  the  potential  of  a  distribution  of  the 
required  type  of  doublets  can  be  obtained  from  a  distribution  of  sources 
located  in  the  plane  of  the  wing. 

The  potential  at  (x,y,z)  due  to  sources  located  at  points  (i,q,0) 
in  region  r  (illustrated  in  fig.  2)  of  the  xy-plane  which  satisfies 
the  differential  equation  (l)  may  be  written  as 


0s 


iOD(x-i) 


COS 


“  R 
M 


R 


d|  drj 


(6) 


where 


& 

it 

and 


Mjl>  M^cjd 


R  =  ^(x  -  |)2  -  (32(y  -  q)2  -  02z2 


Expanding  the  integrand  of  equation  (6)  into  a  power  series  in  to, 
collecting  terms  with  respect  to  |,  and  differentiating  the  resulting 
integral  with  respect  to  z  gives  the  expanded  form  of  the  potential 
of  a  uniform  distribution  of  doublets,  namely 


0D 


Va 

Bz 


2m-3 


+ 


^(all  I  +  a12R  +  •  •  •  +  V2”1’3 

|n(ani  ^  +  an2^  +  •  .  .  +  a-j^R  +  .  .  .j 


d|  dq 


(7) 
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where 


ft 


anl  = 


~  (icD)  ^Ixdx 


ni 


and 


anm  ” 


(■I)”1-1  fcf*-2 

(2m  -  2) 1 \M /  &nl 


(8) 


(9) 


For  convenience  in  the  succeeding  discussion  and  analysis  equa¬ 
tion  (7)  may  he  written  in  the  following  form: 


An  interesting  and  significant  property  of  equation  (10)  is  that 
the  coefficient  of  each  power  of  I  satisfies  the  differential  equa¬ 
tion  (l)  and  has  the  form  of  a  source  potential  with  strength  propor¬ 
tional  to  (“f)ne-la^  •  This  property  may  he  shown  by  writing  the  coef¬ 
ficient  of  |n  as  follows: 


fnl 

R 


1 


(-ir-1  fs 

(2m  -  2)1\M 


j2m-2 


* 
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A  more  general  solution  to  equation  (l)  may  thus  be  obtained  by  intro¬ 
ducing  properly  chosen  weight  or  distribution  functions  (denoted  by 
^Dn(^'n))  into  the  coefficients  of  £n  in  equation  (10).  Let  this 
solution  be  denoted  by  0^)  then  it  can  be  written  as 


* 

* 


Examination  of  equation  (12)  shows  that,  at  the  surface  z  =  0,  the 
potential  0-^  is  determined  by  the  first  integral  expression  but  that 


both  integral  expressions  may  give  rise  to  normal  velocity.  In  succeeding 
steps  in  this  analysis  it  is  shown  that  the  distribution  functions  Dn(t,T}) 


in  equation  (12)  can  be  determined  so  that  the  first  integral  expression 
taken  alone  will  exactly  satisfy  the  boundary  condition  of  tangential 

flow  for  0aj  that  is  fr— — ^  =  Va.  Also,  any  additional  normal  veloc¬ 


ity  that  arises  from  the  second  integral  expression  can  be  canceled,  to 
the  required  order,  by  consideration  of  additional  doublet  solutions  to 
equation  (l) .  The  problem  of  satisfying  the  boundary- value  problem  for 
the  velocity  potential  0a  may  thus  be  reduced  to  that  of  determining 


the  appropriate  distribution  functions  and  additional  solutions  to 
equation  (l). 


In  order  to  show  that  the  first  integral  expression  in  equation  (12) 
can  be  made  to  satisfy  the  boundary  condition  for'  0a,  the  coeffi¬ 
cient  Vaa^  appearing  in  this  equation  is  first  considered.  If  the 
analytical  expression  for  the  coefficient  a^  (equation  (8))  is  multi¬ 
plied  by  xn  and  summed  with  respect  to  n,  the  result  is  identically 
Va.  This  result  may  be  shown  as  follows: 
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It  Is  significant  that  this  identity  holds  if  only  terms  in  55  to  any 
given  power  are  considered.  For  example,  retaining  only  terms  including 
a)  to  the  third  power  gives 


Va 


f  —2  2  ^  ^ 

n  arx^  lapx^ 

1  -  i osx. - 


+  ^  +  X(1S  +  ^ .  a^L\ . 


2X& 


its3 


2  2 


x3  £ 

6 


—  Va 


(14) 


Next  consider  the  normal  velocity  at  z  =  0  associated  vith  0-^,  namely 


W1  = 


,bz  , 


z=0 


Va  lim 
it  z — > 


oSlavl[[  Dn(t,n)  ^  d|  dTi  + 

n=0  JrJ 


y  (*r2anl  r  r 
^(2m-2)lW  JJ 


Dn(l,T1)|n+1R2m-3  dg  dn 

(15) 


Examination  of  this  equation  and  equation  (13)  reveals  that  the  first 
integral  expression  on  the  right  of  equation  (15)  yields  Va  exactly, 
provided  the  distribution  functions  are  determined  so  that  the  following 
integral  equation  is  satisfied: 


Va  lim 
n  z  ?0 


|n+l 

R 


d|  dr)  =  Vax11 


(16) 


The  kernel  of  this  integral  equation  has  the  form  of  a  steady-state 
doublet  potential.  The  problem  of  determining  the  distribution  functions 
for  this  case  is  therefore  analogous  to  determining  distribution  functions 
for  certain  steady-state  problems.  The  distribution  functions  for  steady- 
state  problems,  at  least  for  those  involving  conical  flow,  can  be  deter¬ 
mined  by  a  straightforward  process,  the  main  details  of  which  are  given 
in  the  appendix.  In  this  appendix  a  method  of  solving  equation  (l6)  for 
a  triangular  wing  is  derived  and  the  distribution  functions  required  to 
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derive  the  velocity  potential  for  this  -wing  to  the  third  power  of  <E 
axe  given.  It  is  to  be  noted  that  the  method  derived  for  solving  equa¬ 
tion  (l6)  for  a  triangular  wing  may  be  generalized  to  apply  to  various 
plan  forms  and  to  problems  of  satisfying  the  boundary  conditions  for 
various  velocity  distributions. 

From  this  point  on,  the  analysis  is  restricted  to  the  derivation  of 
terms  of  the  expanded  potential  involving  cn  to  the  third  power.  The 
method  for  deriving  these  first  few  terms  is  quite  general  and  can  be 
used  to  obtain  as  many  additional  terms  of  the  expanded  potential  as  may 
be  desired.  As  previously  pointed  out,  however,  terms  of  the  potential 
after  the  first  few  become  very  unwieldy. 

If  the_appropriate  distribution  functions  are  known  for  terms 
involving  <a  to  the  third  power,  equation  (15)  may  be  written  as  follows: 


Wj  =  Va  +  Vq 


(17) 


where 


w , 


Va  lim  dg 
jt  Z  >0  ^z2 


ixoxJIDqR  + 


d!  dll 


(18) 


is  the  additional  vertical  velocity  arisir^  from  the  second  integral 
expression  in  equation  (12)  involving  m  to  the  third  power.  In  order 
to  maintain  the  boundary  condition  for  0a,  this  additional  velocity  w2 
must  be  canceled.  As  previously  pointed  out  this  canceling,  to  a 
required  order,  can  be  achieved  by  considering  other  doublet  solutions 
to  equation  (l).  For  this  particular  case  consider  a  relation  ^ 

similar  to  (equation  (12))  having  the  following  form: 


T  _  Va  to2  \ 

1  it  dz  2M2  ^0 


nl 


l|n(l>Tl)|n+^  d£  dTl  + 


y~  ( ~i)m~1 

^2  (2“  -  2)  I  \M/ 


nl 


f  f  Dn(  l,T])  |n+%2m“3  d6  dT, 


(19) 
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The  vertical  velocity  w^ ,  to  the  third  power  of  co,  arising  from  this 
expression  is 


wi 


Va  lim  a2  /a5\g 
2*  z - *0  ^z2\M/ 


(1  -  icux)D0  ^—  +  in>Dp  — 
K  R 


drj 


(20) 


In  this  equation  the  distribution  functions  Dq  and  Dp  can  be  deter¬ 
mined, _as  discussed  subsequently,  by  the  method  given  in  the  appendix  so 
that  wp  is  identically  equal  in  value  but  opposite  in  sign  to  w2. 

When  these  functions  are  determined,  the  boundary  condition  for  0a  is 
satisfied  by  wp  +  vp  =  Va,  which  implies  that  the  potential  0  to  the 
third  power  of  co  is  given  by  the  sum 


0a  =  h  +  01  (21) 

to  this  power  of  a>. 

Expressions  for  Dq  and  Dp  are  given  together  with  expressions 
for  other  distribution  functions  D0,  Dp,  D2,  and  D^  in  the  appendix. 
In  regard  to  the  determination  of  Dq  and  Dp,  as  well  as  other  orders 
for  Dn  in  more  extended  treatments,  it  may  appear  necessary,  in  order 

to  formulate  integral  equation#  for  these  functions,  to  perform  the 
generally  unwieldy  integrations  of  the  type 


2^0  f  f  «  *i  (22) 

r 

appearing  in  equation  (15) •  In  general,  however,  the  information  neces¬ 
sary  for  the  determination  of  the  functions  Dn  can  be  obtained,  as  is 
done  in  the  derivation  of  the  functions  Dn  in  the  appendix,  by  exam¬ 
ining  the  values  of  these  integrals  and  their  derivatives  with  regard  to 

the  parameter  6  =  ^  at  some  particular  value  of  G. 

A 

Returning  to  equation  (21)  and  introducing  into  this  equation  the 
expression  for  0p  (equation  (12))  and  the  expression  for  0p 
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(equation  (19)),  each  to  the  third  power  of  o>,  gives  for  the 
potential  0a  to  this  power  of  (5  the  following  result: 


0a-Va2- 


ido  f1  -  ^  -  f *2 + ¥  x3)l 


S2  iai^ 
2M2  2M2 


x  R  + 


t2r>  /.—  ,  —2  ia)3  2\l  io)3  „  ,o_  /aS2  1cd3  \l 

i  Dn  ia)  +  <Xr"X. - X - 5-  R  -  i^Dp - X  — 

1  l  2  /R  2M2  2\2  3  /R 


An  ioA  1.3—  /a^  ia)3  \i  4—  icjp3  1 


6  R 


2M^  2M^  /R 


2M^  R 


Since  it  may  be  shown,  as  in  reference  7,  for  example,  that 


Va  lim  d_ 
jt  z - >0  ct 


Dnln+1  d£  drj  =  VaxnDn(x,y) 


equation  (23)  reduces,  at  z  =  0,  to 


0a  =  Va  1  -  idk  -  ^  X2 


+  x3jDQ(x,y)  + 


x  ia>  + 


^  -  *2(f  -  x)i>2<*»y)  - 


\2  4^-3  \_ 


X3  —•  D3(x,y)  +  x2  j  *—  q  -  ^  x)D0(x,y)  +  x3  D-^y)  (24) 

6  J  \2M2  2M2  /  2M2 
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which,  after  the  expressions  for  the  distribution  functions  given  in  the 
appendix  are  substituted  and  the  terms  are  regrouped,  may  be  written  in 
the  following  simple  form: 


0a  =  Va\/c2x2  -  y2  Aq  +  (A^  -  AQ)iasx  -  cr^a^x2  -  + 


03ia>3x3  +  0^iu>3p2y2x 


(25) 


The  quantities  A q,  Ap,  and  a-^,  0g,  0^,  and  0^  in  this  equation 
are  part  of  a  group  of  quantities  Aj_,  Xj,  ak,  which  are  defined  in 

the  appendix;  these  quantities  are  suitable  for  writing  the  potential 
and  expressions  subsequently  derived  for  forces  and  moments  in  simple  form. 
The  quantities  A±  and  Xj  are  functions  only  of  the  product  0C  (ratio 

of  tangent  of  the  half  apex  angle  of  the  triangle  to  the  Mach  angle)  and 
are  shown  plotted  in  figures  3  and  4,  respectively.  The  quantities  ok 

are  functions  of  Mach  number  M  and  the  product  PC.  These  quantities 
may  be  evaluated  for  particular  values  of  M  from  the  plots  of  Xj  in 

figure  4. 


The  quantities  A q  and  A^  are  the  same,  as  should  be  expected, 
as  the  parameters  associated,  respectively,  with  constant  angle  of  attack 
and  constant  pitching  of  triangular  wings  of  references  8  and  9> 


Expressions  for  the  potentials  0^  and  0£  can  be  obtained  by  the 

method  discussed  for  obtaining  0a,  or  they  can  be  obtained  to  the  third 

order  of  cn  by  comparison  and  synthesis  from  equation  (23).  After 
simplification  these  expressions  are 


0a  =  a-  \fc2x2  -  y^AjX  -  2X2itox2  +  A3inp2y2  -  a^x3  +  06d32p2y2x) - 


Xq\Zc2x2~ T"y2  (a o  -  2Xj_ia»c  -  a^a)2x2  -  agn2p2y2j 


(26) 


0k  =  h 


\jc-x 


2„2 


-  y* 


^Aq  -  2Xjinx  -  a1cu2x2  -  0g a>2p2y2j 


(27) 
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At  (3C  =  1  or  C  =  j-*,  which  is  the  condition  at  which  the  Mach  lines 

from  the  apex  of  the  triangle  coincide  with  the  leading  edges  of  the 
triangle,  equations  (23)>  (26),  and  (27)  reduce,  respectively,  to 


0a  -  p/*2  -  eVfl  -  P  -  ^  5^2  + 

P*  A  3  90M2  90M2 


St±l  1553x3  +  ^2^ 


210M‘ 


210M2 


(28) 


/ 


S  ,  .  §!£  2  215  „2  2  _  gUu£  +  7)  3  . 

,3  *5  *5  315M2 


^(al2  :  7>  -  xii  -  is*  .  Bd  +  2  ^  .  ^2  5^2 


315M* 


90M2  90M‘ 


>t  -  gv*2  -  -  X  -  ^  (30) 


90M 


For  M  =  1  or  for  values  of  the  product  (3C  such  that  p2C2  «  1, 
equations  (23),  (26),  and  (27)  reduce,  respectively,  to 


0a  =  Va/c2x2  -  y2 


(31) 


0a  =  Wc2^  -  y2^x  -  Xq) 


(32) 
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» 

(33) 


These  expressions  axe  the  counterpart  of  the  potential  for  vanishingly 
small  triangular  wings  in  steady  flow  (reference  10 )  and  may  thus  be 
regarded  as  the  potential  of  a  small  triangle  oscillating  in  pitch  and 
vertical  translation  in  either  subsonic  or  supersonic  flow. 


Forces  and  Moments 

As  pointed  out  in  the  introduction,  the  velocity  potential  for  the 
triangular  wing  can  be  used  to  calculate  the  aerodynamic  forces  and 
moments  for  other  plan  forms  that  can  be  formed  from  the  triangular  wing 
by  cutting  the  trailing  edges  so  that  they  lie  ahead  of  the  Mach  cones 
emanating  from  their  foremost  points.  Sketches  of  different  plan  forms 
thus  obtained  are  shown  in  figure  5- 

The  force  and  moment  coefficients  desirable  for  most  flutter  cal¬ 
culations  are  those  that  yield  the  spanwise  variation  in  these  quantities 
or  chordwise  force  and  moment  coefficients.  These  coefficients  are 
obtained  by  integrating  the  pressure  difference  along  any  chord  for  the 
forces  and  the  pressure  difference  multiplied  by  a  moment  arm  for  the 
moments.  A  convenient  procedure  in  deriving  these  quantities  is  to 

introduce  the  reduced- -frequency  parameter  =  k  and  to  employ  the 

the  variables  x,  y,  and  xq  in  a  new  sense  as  nondimens ional  quan¬ 
tities  obtained  by  dividing  the  old  variables  by  the  maximum  chord  2b 
of  the  wing. 

The  pressure  difference  between  the  upper  and  lower  surfaces  of  the 
wing  is 


Ap  =  -2p 


'jL  M 

,2b  dx 


(34) 


The  force,  positive  downward,  at  any  section  of  any  of  the  plan 
forms  shown  in  figure  5  may  be  expressed  as 
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where  has  the  following  values  with  respect  to  the  different  con¬ 

figurations  shown  in  figure  5: 


In  plan  form  A 

H* 

li 

H 

(36a) 

In  plan  form  B 

X-.  =1-2. 

±  m 

(36b) 

In  plan  form  C 

x,  =  1  +  2L 

1  m 

(36c) 

In  plan  form  D 

xn  =  1  +  E 
±  m 

for  0  £  y  <  .  1) 

s 

X1  ~  2bC 

for  m{^c  - x)  <  y g  -k 

> 

(36d) 

After  the  expression  for  0,  given  by  the  sum  of  equations  (25),  (26), 
and  (27),  is  substituted  into  equation  (35)  and  the  integration  is  per¬ 
formed,  the  results  may  be  reduced  to  the  form 


P  =  -i^pbV^e^ 


—  (Lp  +  il^)  +  00(1,3  +  IL4) 


(37) 


where 


Ll  = 


P‘ 


y2/,  2  2  \  Aoy2  ^.1  Cxi 

—  (4-MrX.-^  -  P&Q)  +  — —  cosh  1  — 


(38a) 
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L2  =  \fc2Xl2  -  y2 


T  •  -  2p2^)  ■  |^(3m2c2<,2  +  2^\ ) 


(38b) 


Lo  =  \Zc2Xl2  -  y2-  —  +  ^-(6p2M2X2  -  pVl  -  SM^O!  +  2M2P2X,1)  + 

k2  3P4  ' 


3^2^2A1  "  302M2C2A3  -  3MAC2ct2  -  2M2XxJ 


h  -  SxqLx  (38c) 


^  =  ^2xx2  -  y2 


i-(2p2Ax  +  P2Aq  -  4M2Xx)  +  —  ^  (2^2 


pck 


M‘ 


'p2a1  +  4m^c3  -  4p2M2a5j  -  M  ]^gy  ^4p^2Ag  +  2p2X2  -  M2^  + 


P^C 


4m2P2C2ct2  -  8MAC2ai+  -  8p2C2M2agj 


-  cosh’ 


.1  Cxi 


v 


kC 


y1|‘(4p2C2A3  -  2P2X2  +  M2ax  +  4M2p2C2a2) 


-  2XQL2  (38a) 


In  a  similar  manner,  the  moment  (positive  leading  edge  up)  about 
the  axis  x  =  xq  is 


rx± 

%  =  -4b2  /  (x  -  Xq)  Ap  dx 
Jy/C 


=  -4pV2k2b2 


h0 

—  (Mx  -  iMg)  +  o0(M3  +  mk) 


(39) 
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where 


Mx  =  {c2x±2  -  y2 


-  f>\)  *  +  *\) 


-  2XqL^  (40a) 


Mg  =  -  y2 


Aq  2M2kXn  3 


“^“1  /  V 

xi  + - -r(£P\  -  3M2a1) 


M2kx^2 


p\j 


^p2^  +  M2a1  +  4M2p2C2cr2) 


V; 

Ck 


^2p2Xx  +  M2a1  +  4M2C2a1jj 


cosh 


_i  Cxi 


2xqL2  ( 40b) 


m3  =  \/c2x^2  -  y2 


2 

T~2  (p^A]_  -  4M2p1*'C2A3  -  2p2M2X1  +  2M2P2X2  -  -  4Mi|'p2C2CT2) 

P  c 


‘ioa  ^kh  . 


k2  fr 


p2^  -  6m2P2X2  +  3M^ax) 


“P“  +  ^3  (p4ai  "  4M2p^C2A3  -  2M2p2X1  +  2M2p2k2  -  M4a1  - 


4M^P2C2cr2j 


,  _1  Cxl 
cosh  x  - 

y 


-  2X0^  +  L3  +  2xqL^ 


(40c) 
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M4  =  '/c2x12  -  y 


4xx2 


3P2k 


-  ItM^)  -  ^2j(02Ao  -  S%  + 


2M‘ 


X-J  +  l6M  ^xi.  ^2p\2  -  32M2a1  +  -  4M2p2a5j 


■—  — S- — —  fep^C2A-j  +  23^Xo  -  M2P2a1  +  5M23^2a2  -  M^cr 
153 6C 2  V 


10M^32C2a^  +  M232a5  -  10M23^C2ag)  +  ^|^-{5P6C2A 

4p\2  +  2M232^i  +  5M23^C2a2  +  2M^a3  +  5Mlf32C2a4  - 
2M2p2a^  +  5M2p1|C2ag j  -  2xq^M2  +  L4  +  2xQL2j 


(40d) 


In  equation  (37)^  tor  example,  the  quantity  (L^_  +  il^)  is  the  lift-force 

coefficient  associated  with  vertical  motion  of  the  wing.  The  real  part 
Lp  is  in  phase  with  the  vertical  position  of  the  wing  and  the  complex 
part  L2  is  90°  out  of  phase  with  this  position.  Similar  definitions 
apply  to  the  lift  coefficient  (L^  +  iL^)  associated  with  pitching  motion 
and  to  the  moment  coefficients  (M*j  +  iN^)  and  (M^  +  iMj^) .  The  complex 

or  out-of -phase  terms  determine  the  aerodynamic  damping  associated  with 
different  wing  motions. 

Although  the  expressions  for  the  components  of  lift-force  and  moment 
coefficients  in  equations  (38)  and  (40),  respectively,  are  lengthy,  they 
may  he  quite  easily  evaluated  with  the  aid  of  the  graphs  in  figures  3 
and  4. 
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DISCUSSION 
Sample  Calculations 


In  order  to  give  some  indication  as  to  the  general  nature  of  the 
spanwise  distribution  of  the  different  components  of  lift  and  moment 
coefficients,  equations  (38)  and  (40)  have  been  evaluated  at  different 
spanvri.se  positions  y  for  plan  form  A  (fig-  5(a))  for  ‘the  following  set 

of  conditions:  f3C  =  0.5,  C  =  1.0,  xq  =  0.6,  M  =  ^5/4 >  and  k  =  0.1. 

These  sample  results  are  plotted  as  functions  of  spanwise  position  in 
figure  6.  The  spanwise  variations  of  the  different  components  of  lift 
force  are  shown  in  figure  6(a)  and  the  corresponding  variations  of 
moment  coefficients  in  figure  6(b) . 

In  figure  6  note  that,  for  the  particular  set  of  conditions  for  this 
example,  the  maximum  values  of  the  components  of  moment  coefficients  M^_, 

M2,  and  M3  are  positive  and  act  near  the  tips  of  the  wing,  whereas  the 
maximum  value  of  the  component  Mj^_  is  negative  but  also  acts  near  the 

tips.  It  may  also  be  noted  that  the  integrated  (in  spanwise  direction) 
values  of  the  components  of  moment  coefficient,  or  components  of  total 
moment  coefficient,  would  in  each  case  have  the  same  sign  as  the  maximum 
value  of  the  corresponding  component  of  section  moment  coefficient.  This 
result  is  not  necessarily  true  in  general,  because  changing  some  of  the 
parameters  involved  in  the  evaluation  of  the  spanwise  distribution  of 
some  components  of  both  force  and  moment  coefficients  may  change  the 
distributions  significantly  from  those  shown  in  figure  6. 

The  fact  that  the  total  component  of  moment  coefficient  M^  is 

negative  in  the  example  just  discussed  shows  that,  for  the  conditions  of 
the  example,  this  term  would  not  contribute  to  the  aerodynamic  damping 
but,  on  the  contrary,  would  act  as  a  source  of  energy  for  the  oscillating 
system.  This  circumstance  is  significant  since  it  leads  to  the  possi¬ 
bility  of  the  single-degree-of -freedom  torsional  instability  discussed 
in  the  following  paragraphs. 


Undamped  Torsional  Oscillations 

The  wing  plan  forms  discussed  herein,  like  two-dimensional  and 
rectangular  wings,  exhibit  the  possibility  of  undamped  torsional  oscil¬ 
lations  for  certain  ranges  of  Mach  number  M  and  location  of  axis  of 
rotation  xq.  This  fact  is  borne  out,  as  indicated  in  the  preceding 
paragraph,  by  considering  the  integrated  (spanwise)  value  of  the  com¬ 
ponent  of  damping  moment  Mj^  associated  with  pitching  or  torsional 
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motions.  The  main  results  of  this  phenomenon  can  he  obtained  by  con¬ 
sidering  very  slow  oscillations  so  that  only  terms  in  equation  (40d) 
for  M]^  involving  the  reduced  frequency  k  to  the  order  l/k  need  be 

retained.  In  this  case, 


+ 


8(2M2  -  1)  2 

3P2C2  17 


+ 


(2M2  -  l)A1)x1  +  (BSxqAq 


(to) 


For  the  triangular  plan  form  (fig.  5(a))  the  integrated  value  or  total 
component  of  moment  coefficient  is 


m4  = 


JtC 


202k 


-(2M2  +  1)Aq  +  3(2M2  -  1)AX  +  ^[Aq  -  (2M2  -  ljA-H  + 


toQ2(M2  -  1)Aq 


(to) 


In  general,  the  condition  of  torsional  stability  or  instability 
depends  on  the  sign  and  magnitude  of  in  equation  (42).  Positive 

values  of  indicate  stable  conditions  and  negative  values  indicate 

the  possibility  of  torsional  instability.  Between  the  stable  and  unstable 
conditions,  that  is  when  vanishes,  a  borderline  state  of  unstable 

equilibrium  separating  damped  and  undamped  torsional  oscillations  exists. 

The  ranges  of  valuesjof  Mach  number  M  and  location  of  axis  of 
rotation  Xq  for  which  vanishes  for  some  selected  values 

of  C  =  tan  6  are  shown  plotted  in  figure  J.  The  regions  inside  the 
curve  in  this  figure  indicate  instability.  The  dashed  curve,  on  which 
some  of  the  solid  curves  terminate,  represents  the  locus  of  values  of  M 
and  xq  for  which  vanishes  for  the  whole  class  of  triangular  wings 

with  supersonic  leading  edges,  that  is,  for  triangular  wings  where 

PC  >  1.  It  will  be  noted  that  vanishes  for  values  of  Xq  ahead  of 

the  root  2--  chord  position.  It  will  also  be  noted  that,  as  the  vertex 
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angle  e  =  tan“^  C  decreases  to  30o>  the  range  of  values  of  Mach  number 
for  which  vanishes  decreases  sharply. 

In  conclusion,  investigation  of  equation  (42)  shows  that,  for  a 
given  value  of  the  reduced  frequency  k,  Mach  number  M,  and  location  of 
the  axis  of  rotation  Xq,  the  magnitude  of  the  damping  coefficient 

decreases  as  C  decreases  and,  consequently,  torsional  instability  is 
less  likely  to  occur  with  slender  triangles  than  with  wider  triangles. 


Langley  Aeronautical  Laboratory 

National  Advisory  Committee  for  Aeronautics 
Langley  Field,  Va.,  June  19,  1951 
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APPENDIX 


METHOD  FOR  DETERMINING  DISTRIBUTION  FUNCTIONS 


In  this  appendix  a  method  of  obtaining  the  distribution  functions 
is  developed  in  detail.  Expressions  for  the  distribution  functions 
required  to  derive  the  velocity  potential  to  the  third  power  of  the 
frequency  of  oscillations  and  a  set  of  functions  useful  in  writing  the 
expression  for  the  potential  in  concise  form  are  given. 

As  indicated  in  references  8  and  9y  where  the  triangular  wing  is 
treated  for  constant  angle  of  attack  and  for  constant  rolling  and  pitching 
motions,  a  convenient  form  of  the  velocity  potential  corresponding  to  a 
distribution  of  vertical  velocity  proportional  to  xn  is 

K  =  “  da 

or  dz 


L 


tn+2  a  i 

i  d| 


0  \fU  -  I)2  -  P2(y  -  a!)2  -  P2z2 


(Al) 


where  a  is  the  slope  of  a  ray  passing  through  the  vertex  of  the  wing 
(*1  =  01);  !i  is  the  least  value  of  |  that  causes  the  denominator  R 

in  the /integrand  to  vanish,  and  Fn(a)  =  j  Dn(£,r))  is  the  distribution 
function  that  is  to  he  determined  so  that 


(wn) 


z=0 


(— )  =  - 

'dz  4=0 


(A2) 


In  equation  (Al),  the  integrand  is  noted  to  be  singular  at  the 
limit  £  =  A  form  of  the  integrand  which  avoids  this  difficulty  is 

obtained  by  making  the  following  change  of  variables: 


q 


X 

1  - 


(1  -  P2 aQ) 


N  = 


px 


1  -  p2a2 


'(«  - 


„)2  +  %  (1  -  (S^) 


0  =  £ 
X 


T  =  cosh' 


■1  Q  -  g 
N 


(A3) 


where  t  is  the  new  variable  of  integration.  With  these  substitutions  equation  (Al)  becomes 
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passing  to  the  limit  gives  for  equation  (A5) 
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w  ,  P2(n  +  2)  rC  ?n(»)  3°  rfash-1  A 

n  n  '“  -C  (l  -  p2a2p/2 

)  ,  *1 
cosh  T  dT 

=  Xn  (A 6) 

where  Nq  is  the  value  of  N  at  z  =  0.  The  requirement  that  the  normal 

velocity  be  proportional  to  xn  and  independent  of  y  implies  that  all 
derivatives  of  equation  (a6)  with  respect  to  y  or  0  must  vanish,  or 
that  Fn(a)  must  be  so  determined  that  the  final  value  of  the  integral 
in  this  equation  be  independent  of  6 .  The  requirement  that  all  deriv¬ 
atives  of  equation  (A6)  with  respect  to  6  vanish  leads  to  other  equa¬ 
tions  for  Fn  and,  after  n  +  1  such  differentiations,  the  equations 

acquire  forms  for  which  solutions  are  known.  The  value  of  Fn  can  then 
be  determined  from  these  known  solutions  by  evaluating  each  of  the 
aKvn 

derivatives  — —  at  any  arbitrary  value  of  6  in  the  range  -C  <  Q  <  C 

ae* 

as  follows: 

The  Kth  derivative  of  vn  with  respect  to  6  ( equation  (a6) )  gives 
the  following  integral  equation  for  Fn: 


d0K 


P2(n  +  2) 

n 


Fn(o)  da  f" 

(T77vpJ0 


Nn  cosh  T^n+^  — -  <t(  T )  dT 

u  >  N0 


0 


(A7) 
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where 


0(t)  = 


(n  +  l) !  /  Pa  +  cosh  t 


iK 


(n  +  1  -  K)! \N0  cosh  t  -  q 


K(n  +  l) !  /  Pa  +  cosh  t  \^-l 

(n  +  2  -  K) !  \Nq  cosh  t  -  ql  ® 


K(K  -  l)(n  +  1):/  Pa  +  cosh  t  \K~2tt  -2 
( n  +  3  -  K) !  \N0  cosh  t  -  q) 


V  + 


.  +  (-l)KK!Nn"K 


When  K  =  n  +  1,  the  expression  for  4>(t)  may  be  recognized  as  being  a 
binomial  expression,  namely 


$(t)  as  (n  +  l)  I 


P a  +  cosh  t 
Nq  cosh  t  -  q 


i_\n+1 

No/ 


(n  +  1) !(paN0  +  q)n+1 
(-No)n+1(q  -  Nq  cosh  t)11^ 


_ (n  +  l)  !xn+1 _ 

(-No)n+1(q  -  Nq  cosh  T)n+1 


(A8) 


Thus  the  integral  equation  corresponding  to  the  n  +  1  derivative  of 
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This  function  alBO  satisfies  equation  (A9)  for  any  value  of  n  and 

d-K-W-n 

satisfies  equation  (A7)  for  all  derivatives  — beyond  the  nth 

d@K 

derivative.  A  more  general  form  of  solution  to  equation  (A10)  may  be 
shown  to  be 


m 


(A12) 


where  m  is  an  integer  and  m  ^  n.  The  validity  of  this  solution 
follows,  since 


r  i  m 

0®  =  -  (6  -  ff) 


m!em-r 


(m  -  r) !r! 


(0  -  o)r 


(A13) 


and 


Wc2  -  g2  _  miem~r  fc2  -  a2 

(0  -  cr)n+3  ~  rR)  (m  -  r)lr!  (6  -  a)n’r+3 


(Al4) 


Each  term  of  the  summation  in  equation  (Alt)  is  found  to  have  the  form 
of  the  integrand  in  equation  (A10),  from  which  the  function  (A12)  is 
concluded  to  satisfy  equation  (A10)  and  equation  (A7)  for  all  values  of 
K  >  n. 

From  the  foregoing  discussion  and  consideration  the  distribution 
function  Fn(cr)  may  be  uniquely  determined  in  terms  of  expression  (A12) 
as  follows: 

Consider  the  expression 


Fn(ff)  “  A0f0  +  Alfl  +  •  •  •  +  Anfn  (A15) 


where  the  coefficients  A q,  A-^,  .  ,  .  An  are  constants  that  are  to 
be  determined. 

Each  term  on  the  right  of  equation  (A15)  is  noted  to  satisfy  equa¬ 
tion  (A7)  for  all  derivatives  beyond  the  nth  and  a  total  of 
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n  +  1  parameters  is  to  be  determined.  If  expression  (A15)  is  intro- 
duced  into  equations  (A6)  and  (A7)  and  the  indicated  integrations  in 
equation  (a6)  and  in  equation  (A7)  for  K  =  1,  2,  .  .  .  n  are  performed, 
n  +  1  linear  equations  are  obtained  in  Aq,  A^,  .  .  .  An  from  which 

these  constants  may  be  determined. 

The  integrations  with  respect  to  a  in  equations  (A 6)  and  (A7)  are 
in  general  difficult  and  tedious  to  perform;  however,  as  previously 
pointed  out  it  is  only  necessary  to  perform  the  integrations  for  some 
particular  value  of  0  in  the  range  -C  <0  <  C .  The  integrals  have 
their  simplest  form  when  0  =  0  and  the  integrations  can  be  made  for 
this  value  of  0  by  reductions  and  use  of  formulas  in  reference  11. 

The  distribution  functions  Fq(cx),  F^(cr),  FgCcr),  F^cr),  Fq(cj), 
and  F^(a)  and  functions  Bn(|,r])  calculated  by  the  foregoing  method  are 
as  follows: 
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a2  = 


a3  = 


a4  = 


A 6  = 

a?  = 

A 6  = 

^9 


2d3 

2d]_dg 

+  p^d^ 

2d2 

2d]_d3 

+  p2C2d22 

6d6 

6d^d6 

+  P2C2d52 

6d5 

6d4d6 

+  P2C2d^2 

2d3d7 

-  P2C2d2dg 

2d^d-^ 

+  p2C2d22 

2(dld8 

;  +  d2d7) 

2d^3 

+  P2C2d22 

6d6d9 

-  P2C2d2d^ 

6dl^dg 

+  p2c2d52 

6( dgd^ 

+  d^d^) 

6d4d6 

+  P2C2d52 

The  quantities  dn  are  defined  as 


5(P^C2  -  3P1<C4)F'  +  (2  -  lOp^2  +  gp^E1 
2(1  -  p^2)2 
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2p2C2F'  -  (1  +  p2c2)E’ 
(1  -  P^2)2 


-  ^pV^F'  -  (6  -  lOP^2  +  2p*tC^)Et 
(1  -  p^2)2 


(27p2c2  -  sipk:1*-  +  12P^c^)f’  +  (6  -  ^p2^  +  65$^_  -  24p^c^)E' 

6(1  -  p2c2)3 


(9p2c2  -  p^)F'  -  (3  +  7P2C2  -  2p^)E' 

(1  -  P2C2)3 


(3p2c2  _  pp1^  _  2P^c6)F'  -  (6  -  15p2c2  +  5Pi,C^  -  I+p6c6)E' 

(1  -  p^2)3 


P2C2(E'  -  F1) 

2(1  -  p^2) 


p2c2j"  _  e' 

1  -  p2c2 


P2C2[(P2C2  -  3)F'  +  (Ij-  -  2P2C2)E'_ 

6(1  -  p^2)2 


5 
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Combinations  of  these  functions  useful  for  writing  the  potential 
in  concise  form  are : 


*0  -  *1  .  X, 

2  1 

*5  " 

A1  -  a2 

2  -  *2 

a3 

2 

A2  "  A4  _  y 

2  3 

A 7 

3 

Xl  -  X2  =  *5 

a8 

2 

^  -  *3  =  *6 

II 

1 

x6  - 

A3 

*4  +  2  =  *8 

X^  + 

A(^6  , 

2  -  *9 

A0A7  , 

2  =  *10 

A1A8  -  Am 
—£ —  =  K11 

aia9  , 

— 2~  ~  *12 

X9  -  An  =  *13 

*9 

M2 


10 

^“=CT2 


*0.1 

M2 


=  Ot 


■  *12  _ 
M2 


°6 


*10  -  A12  =  *14 


It  will  be  noted  that  the  quantities  A±  (i  =  0,  1,  ...  9)  and 
Aj  (j  =  1,  2,  .  .  .  l4)  are  functions  only  of  the  product  pC .  They 

are  shown  plotted  in  figures  3  and  k,  respectively.  The  quantities 
(k  =  I*  2,  .  .  .  6)  are  functions  of  Mach  number  M  and  the 

product  pC  and  can  he  evaluated  for  particular  values  of  M  from  the 
plots  of  A±  and  Xj. 
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;etch  illustrating  region  of  integration  for  the  velocity 
potential. 
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(a)  j  =  1,  2,  and  3. 
Figure  4.-  Variation  of  the  quantities  X.j 


pc. 


as  functions  of 
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(U) 


j  =  4,  5,  6,  7,  8,  and  9* 
Figure  4. - 


L 


Continued 


4o 
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(c)  j  =  10,  11,  12,  13,  and  l4. 
Figure  4.-  Concluded. 


(t>)  Moment. 
Figure  6.-  Concluded. 
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NACA-Langley  -  9-17-51  -  1000 


Figure  7*-  Curves  shoving  ranges  of  Mach  number  M  and  axis  of 

rotation  x0  for  vhich  the  aerodynamic  torsional  damping  moment 
vanishes  for  some  selected  values  of  angle  e. 


* 


tr¬ 

io 

CM 


1 

53  JS  S 

ll-g 

U  0)  S 
ftg  rt 

'S  a§ 
«£S 
a  ■g* 

Sod 
g  Xh  o 

<D 

0)  >  g 

u  0)  <D 

&3  “ 

*  S.2 

m  *«  ? 
O  ca 

^  HI  h 
+$  >  O 

8  •Eb’g 


g  rt  w  « 

“  w"S  c 
CO  c  ®  .Q> 

•f-i  O  C  *Cl 

.fig  §,£ 

a  3  fi 

■gS  o  o 

?  o.  o  -o 


■g 


9«1 

•3  S  6 

S3a 

■si 

* 

«*S 

SH" 

'S  Tj  ° 

Ql  T-t 
QJ  >  -M 

»*  g 

<33  w 
^  .  0) 
to  £  w 

s-a-p 

•2  g-g 

Stt-g 

O  0)  w 
^  h  O 

g  rt  to  w 
“  wTc  c 

a  §  §2 

xi  ^  §,32 

fuss 


12  8 

!§ 
rt  a 

•H  w 

g  T3 

8  g 

-M  TO 

■gftS 

tai^S 

a  -a* 
^•33 

&>* 

?7  <D  <D 

33  “ 

^  <J) 

rt  S  <2 

^  *c  *§ 
-e  >  o 
§  ’Sb'o 

§  9  •* 

*3  b  o 

MH  TO  *n 


2  ®  «r-s  c 

P  .2  §§.2 


<J  X! 
1  * 


*s 


C  C 
<U  0) 
e  .h 
§  « 

fi 

.  o  o 
a  o  o 


53 


•  +, 

S-o  g 

« s  g 

s 

c  -a 
a>  3 

-4-»  TO 

o,  °> 
^  o 

£>S 
*3  « 
3  g 

a>  x. 

^  o 
a>  a? 
xs  w 


n*2? 

O  c  £ 
^  a>  »h 

14! 

223 

*+h  TO 


E  « 


CQ  W 

~-4-»  •*-» 

W  C5  C 
C  QJ  0) 
O  C  7* 
rj  O  .3 
rt  a«fl 

8  So 

32  o  o 
0,00 


% 


* 


